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I.  ^INTRODUCTION 

This  paper  gives  a  brief  description  of  a  shortcut  method,  the  H-R^ 
method,  for  predicting  transition  in  a  vide  class  of  boundary  layer  flows, 
including  the  effects  of  pressure  gradient,  surface  heat  transfer  and  sue- 
tion.  Here  H  and  R^/are  the  body  shape  factor  and  the  Reynolds  number, 
based  on  distance  x  measured  in  the  direction  of  the  flow,  respectively. 
The  method  is  extremely  simple  to  use  and  a  good  substitute  to  the  well 


known  but  rather  complicated  *.  method.  ( <  -  *' 

K  "  f  4 ' 

II.  TOE  e  FETHOD  OF  FORECASTING  TRANSITION 

The  H-R^  method  has  not  been  correlated  with  test  data  but  rather 

Q 

has  been  justified  in  terms  of  Tollmien-Schlichting  waves  and  e  type  cal- 

Q 

culations.  Data  are  gradually  being  accumulated  that  show  that  the  e 
method  is  the  best  all  around  method  that  now  exists  for  predicting  boun- 

Q 

daxy  layer  transition.  Since  the  H-Rx  method  is  rooted  in  the  e  method, 

we  proceed  to  make  some  comments  about  the  latter. 

Transition,  although  it  may  commence  with  the  amplification  of  Tollmien- 

Schlichting  waves  as  described  by  linear  stability  theory,  is  dominated  in 

its  late  stages  by  three  dimensional  and  non-linear  effects.  Why  then  does 

0 

transition  occur  at  a  disturbance  amplification  ratio  at  of  about  *  e 
[1]  as  computed  from  linear  instability  theory? 

9 

The  e  method  is  rooted  in  the  following  observations:  Liepmann  [2] 

hypothesized  that  at  the  breakdown  to  turbulence,  the  Reynolds  stress 

T  ■  -pu'v'  due  to  the  amplified  fluctuations  becomes  comparable  in  magni- 

*  Consultant  The  Rand  Corporation, Santa  Monica, California;  also  Professor  The 
University  of  California, Los  Angeles. 

•♦Consultant  The  Rand  Corporation, Santa  Monica, California;also  Adjunct  Professor 
The  University  of  California, Los  Angeles.  Associate  Fellow  AIAA. 

***  Adjunct  Professor  The  University  of  California, Los  Angeles. Fellow  AIAA. 


tude  to  the  maximum  mean  laminar  shear  stress  =  Vi(3y/3y).  The  ratio  t/t^ 
is  given  by 

T/T,  -  ^  [kb^Vw].^  (1) 

where  c,  is  the  skin  friction  coefficient,  u'  is  the  disturbance  amplitude 

*L  _ n 

at  the  neutral  point,  b  *  v'/u'  and  k  =  u'v’/uv,  a(x)  is  the  amplification 

factor  with  respect  to  the  neutral  point,  u  and  v  are  velocities  in  the  x  and 

y  direction  respectively  and  prime  indicates  disturbance  values.  Obremski 

et  al.  [3]  observes  that  in  a  low  turbulence  environment  the  critical  Tollmien- 

Schlichting  mode  at  the  beginning  of  ampl if i cation  may  possess  an  amplitude 

u'  of  the  order  of  0.001%  u  ,  or  perhaps  even  less;  u  is  the  edge  velocity, 
n  c  c 

Furthermore,  Klebanoff  et  al.  [4]  found,  in  detailed  examination  of  flat 
plate  measurements,  that  disturbance  growth  via  the  laminar  instability  mecha¬ 
nism  ceases  to  be  valid  when  the  rms  velocity  fluctuation  v'  in  the  boundary 
layer  reaches  (u'/u  )  s  0.015,  but  that  the  first  appearance  of  turbulent 

v 

spots  is  expected  at  about  (u'/u.)  =  0.20.  Assuming  that  these  figures 

represent  disturbance  growth  not  only  on  a  flat  plate  but  aLt>o  on  airfoils 

and  bodies  of  revolution,  all  at  low  free  stream  turbulence  level,  we  find 

that  amplification  in  the  linear  regime  [(u'/u  )  /(u'/u  )  z  0.015/0.00001  * 

c  max  n  c 

7  35 

1S00  sc*  ]  substantially  exceeds  the  amplification  in  the  nonlinear  regime 

[0.20/0.015  c  13  s  e2'**].  Total  amplification  (*  0.20/0.00001)  at  transition 

9  9 

is  seen  to  be  about  20,000  s  e  ’  ;  this  amplification  is  of  the  same  order  as 

9  2 

that  reported  by  Michel  [5],  e  ’  ,  and  computed  by  Smith  et  al.  [6]  and  van 
Ingen  [7]  from  linear  instability  theory.  This  amplification  factor  is  con¬ 
sistent  with  the  hypothesis  of  Liepmann.  If  we  set  t/t^  *  1,  (u^/ug)  = 

0.00001,  cf  s  0.664Rl“1^2  with  =  3xl06  (for  transition  on  a  flat  plate), 

L 

(u'/v')tr  =  0(10)  and  (uv/urvT)tr  *  0(10)  at  transition,  then  eq.  (1)  gives 
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an  amplification  at  transition  *(*tr)  =  1.4x10  ■  e  This  suggests  that 


in  boundary  layer  flows,  where  free-stream  turbulence  is  very  low,  the 
growth  of  Tollmien-Schlichting  waves  controls  the  development  to  turbulent 
flow  and  so  linear  instability  theory  can  be  used  as  a  basis  for  forecast¬ 
ing  transition. 

III.  S01€  PREDICTIONS  FROM  STABILITY  THEORY  AND  THE  e9  KCTHOD 

The  factors  being  considered  in  this  note,  for  either  two-dimensional 
or  axisymmetric  low  speed  flow,  are  (1)  effect  of  pressure  gradient,  (2)  ef¬ 
fect  of  suction,  and  (3)  effect  of  wall  heating  or  cooling.  Figure  1  shows 
the  results  of  calculations  of  the  neutral  stability  point  for  wedge  and 
other  flows  under  a  large  variety  of  conditions.  It  contains  results  for 
an  adiabatic  flat  plate  with  varying  degrees  of  mass  transfer,  including 
the  asymptotic  suction  case.  The  largest  number  of  points  are  for  various 
wedge  flows  in  water  with  and  without  heat.  The  significant  fact  is  that  the 
neutral  stability  point  when  plotted  in  the  fox*  of  the  critical  Reynolds 
number,  Rg*crit»  vs-  H»  the  shape  factor,  this  rather  considerable  variety 
of  flows  forms  a  6-ingZe.  well  defined  curve.  Shape  factor  is  the  immediate 
determinant  of  neutral  stability  rather  than  some  more  remote  measure  such 
as  Hartree's  3  or  Pohlhausen's  A. 

The  present  authors  proceeded  to  determine  if  a  similar  correlation 

g 

existed  at  transition  as  predicted  by  the  e  method.  Suction  was  not  in¬ 
cluded  in  the  study  but  a  considerable  number  of  wedge  flows  with  various 
wall  temperatures  were  studied  for  water.  The  results  are  shown  in  Figure 
2,  which  also  includes  the  data  of  Figure  1.  But  in  Figure  2  the  calcu¬ 
lated  values  are  plotted  in  terms  of  Rx  vs.  H  rather  than  Rg*  vs.  H.  Again, 
except  for  the  highest  temperature  differentials  the  data  form  well  defined 
curves;  the  scatter  is  so  small  that  suction  effects  would  undoubtedly 
fall  along  the  same  curve.  The  correlation  indicates  that  the  neutral  and 
transition  points  are  primarily  a  function  of  H.  Pressure  gradient,  heating 


and  suction  are  only  methods  of  influencing  H.  The  results  of  Figure  2  are 

only  for  two-dimensional  and  similar  flows,  but  they  suggest  a  method  of 

predicting  transition:  for  the  conditions  of  the  problem,  make  a  plot  of  H 

vs.  R^,  from  ordinary  boundary  layer  calculations,  as  one  proceeds  back  along 

the  body.  The  plot  will  always  start  beneath  the  two  loci  of  Figure  2. 

g 

When  the  curve  crosses  the  e  locus,  transition  should  occur.  An  equation 

9 

that  fits  the  e  locus  very  will  is  (for  bodies  of  revolution  R^  is  replaced 
with  Rs  where  s  is  the  distance  measured  along  the  surface  of  the  body) : 

log[Rx(e9)]  -  -40.4557  ♦  64.8066H  -  26.7538H2  ♦  3.3819H3  ,  (2) 

2.1  <  H  <  2.8 

IV.  THE  H-Rx  METHOD  FOR  PREDICTING  TRANSITION  OVER  A  BODY  OF  REVOLUTION 
Because  the  correlations  of  Figures  1  and  2  are  for  two-dimensional 
similar  flows,  as  a  sample  problem  we  chose  a  heated  Reichardt  body  of  revo¬ 
lution  moving  in  67°F,  19.44aC,  water  at  a  velocity  that  gives  um/\>  =  6xl06/ 
ft,  20xl06/m,  Figure  3.  Because  there  is  a  stagnation  point  the  flow  starts 
out  very  stable  (low  H) .  For  the  case  of  no  heating  (AT  =  0°C)  the  neutral 

point  is  crossed  at  R^  z  10 ®  and  transition  occurs  between  R^  =  3  and  4xl06. 

g 

The  circle  is  for  the  full  e  calculation,  so  the  agreement  is  perfect.  With 
10"F  (5.5°C)  of  heating,  the  boundary  layer  is  more  stable  and  now  the  neutral 
point  is  crossed  at  about  Rx  *  1.2xl06  and  transition  is  moved  back  to  about 
Rx  *  9x10^.  Again  the  agreement  with  the  full  e9  method  is  very  good.  At 
a  heating  of  20#F  (11.1°C)  transition  is  moved  still  further  back  but  the 

9 

agreement  with  the  full  e  calculation  is  not  quite  as  good.  But  the  example 
shows  clearly  the  method  of  prediction  and  that  the  correlation  seems  to  work 
even  when  the  flow  is  nonsimilar  and  non- two-dimensional . 

For  a  practical  shape  such  as  in  this  example,  the  H-Rx  trace  meanders. 
For  a  pure  similar  flow  the  trace  would  be  vertical.  If  there  is  a  sudden 


change  in  pressure  gradient,  suction  or  heat  the  trace  would  be  horizontal 
for  some  distance.  Note  that  only  ordinary  but  accurate  boundary  layer  cal¬ 
culations  of  H  vs.  Rx  are  needed  for  transition  prediction,  not  lengthy 
stability  calculations. 

V.  DISCUSSION 

Other  simple  methods  of  predicting  transition  lack  the  generality  of 
the  H-R^  method  because  they  are  correlations  based  on  parameters  that  are 
more  removed  from  the  best  measure  of  stability,  H.  Michel's  method  for 
instance  is  an  excellent  correlation  method,  correlating  Rg vs .  R%  at  transi¬ 
tion  for  two-dimensional  incompressible  flow,  e.g.,  airfoils.  But  when 
bodies  of  revolution,  heat,  or  suction  are  considered  the  correlation  is 
no  longer  applicable.  The  same  kind  of  trouble  applies  to  methods  like 
Granville's  [8],  whibh  correlates  with  a  pressure  gradient  parameter  such 
as  Pohlhausen's  A.  Obviously  this  kind  of  method  fails  if  other  means  than 
pressure  gradients  are  used  to  affect  transition.  Consider  a  fixed  point 
on  a  particular  body.  That  point  has  one  value  of  A  regardless  of  the  degree 
of  suction  or  heating.  But  if  H  is  used  as  the  correlating  parameter,  changes 
at  that  point  in  the  boundary  layer  flow,  as  by  heating  or  suction,  manifest 
themselves  in  changes  in  H.  However,  we  must  stress  that  when  there  is  heat 

transfer  more  than  the  shape  factor  is  involved,  particularly  if  temperature 

q 

differences  are  high.  For  instance,  in  the  full  e  method  using  the  extended 

2  2 

Orr-Sommerfeld  equation  [1],  additional  terms  such  as  d  y/dy  enter  the  sta¬ 
bility  problem  in  addition  to  the  basic  velocity  profile  data  that  determines 
H.  Note  that  the  correlation  is  poor  at  high  heating  rates  where  H  <  2.1, 
Figure  2. 

The  H-R^  method  for  predicting  transition  gives  reasonable  answers  so 
long  as  the  flow  does  not  vary  too  much  from  nearly  similar  flow,  i.e.,  local 


5 


similarity.  Also  it  is  applicable  only  so  long  as  the  effects  of  surface 

rou^mess,  vibrations  and  freestreaa  turbulence  level  are  sufficiently  low, 

q 

just  as  required  for  the  basic  e  method.  The  method  is  also  restricted 

to  heating  rates  where  T  -T^  does  not  exceed  about  23*C. 

The  correlation  has  been  developed  entirely  theoretically  using  the 
9 

e  method  as  if  it  were  exact- for  predicting  transition.  No  attempt  has 

been  made  to  develop  the  correlation  from  experimental  data,  partly  because 

of  the  labor  involved  but  partly  because  good  test  data  are  quite  scarce 

for  unusual  conditions  such  as  heating  in  water.  The  main  fact  that  can 

be  claimed  therefore  is  that  the  H-Rx  method  is  a  very  convenient  substi- 

9 

tute  for  calculating  nearly  the  same  results  as  the  full  e  method.  Also 
of  course  this  saiqtle  case  is  not  the  only  one  that  has  been  studied;  the 
H-Rx  method  has  been  used  extensively  in  other  studies  with  plausible  predic 

9 

tions,  and  a  number  of  other  comparisons  with  the  full  e  method  have  been 

made.  Where  the  flow  differs  widely  from  the  conditions  of  the  correlation, 

9 

a  check  should  be  made  by  the  full  e  method.  But  for  a  large  range  of 
practical  conditions,  the  H-R^  method  seems  accurate,  convenient  and  quite 


general . 
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Critical  Reynolds  number  vs.  H. 

Critical  and  transition  (computed)  Reynolds  number  for  adiabatic 
and  heated  wedge  flows  in  water.  $  is  Hartree's  6 . 

Paths  of  boundary- layer  development  (H  vs.  R^)  and  predicted 
transition  for  a  13:1  Reichardt  body.  Conparison  with  heated 
wedge-flow  predictions.  Circles  denote  point  where  Re^^  is 
reached,  u^v  ■  6xl0^/ft,  20xl0^/m.  s  is  length  measured  along 
the  body  surface. 


